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We study the effect of spin-orbit scattering on the statistics of the conductance of a quantum dot 
for Coulomb blockade peaks and valleys. We find the distribution function of the peak heights for 
strong spin-orbit scattering in the presence and absence of time reversal symmetry. We find that the 
application of a magnetic field suppresses the average peak height, similar to the antilocalizaion in the 
bulk systems. For the valleys, we consider the elastic cotunneling contribution to the conductance 
and calculate its moments at the crossover between ensembles of various symmetries. 
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Electron transport in quantum dots (QD) is affected by 
electron-electron interactions, the profound example be- 
ing Coulomb blockade (CB) l] in QDs connected to leads 
by contacts of small conductance G /ttH = Gq. Due 
to charge quantization, at temperatures smaller than the 
charging energy i?c ^ A (A being the one particle level 
spacing), associated with the addition of one electron to 
the QD, transport through the dot is diminished except 
at discrete values of the gate voltage corresponding to 
charge degeneracy. Thus, the dependence G{Vg) exhibits 
almost equidistant sharp peaks separated by deep min- 
ima (CB valleys). Close to the peaks and at low tem- 
peratures T ^ A, transport occurs by resonant tunnel- 
ing through a single state. The situation is different in 
the valleys where transport is due to virtual transitions 
of electrons via excited states of the QD: many levels 
(TV ~ Ec/A) contribute to tunneling. In both cases, 
the mesoscopic fluctuations of the conductance are de- 
termined by statistics of wavefunctions in the closed dot. 

In disordered or chaotic QDs with the Thouless en- 
ergy i?Th much larger than all relevant energy scales, 
transport phenomena are described by Random Matrix 
Theory (RMT) (see e.g. 0). In this theory, the statis- 
tics of one particle energy levels and eigenfunctions are 
described by universal ensembles sensitive only to the 
underlying symmetries. The statistical distribution of 
the peak heights has been calculated based on RMT for 
the Gaussian orthogonal (GOE) and unitary (GUE) en- 
sembles 13]. The results were extended to the crossover 
between those ensembles 0, using the statistics of wave- 
functions derived in Ref . |a. The average peak height was 
observed to increase with the magnetic field, reminiscent 
of weak localization in bulk systems [^. The statistics 
of conductance in the valleys was studied theoretically 
and experimentally Q for the same crossover. How- 
ever, there is a gap in the theoretical literature regarding 
the effect of spin-orbit scattering on these statistics. The 
effect of breaking of time reversal symmetry on the CB 
peak heights in quantum dots with strong spin-orbit scat- 
tering is expected to be strong, similar to antilocalization 
in bulk systems or open dots [lo| . 

In this Letter, we find the statistics of the peak heights 



for strong spin-orbit scattering in the presence and ab- 
sence of time reversal symmetry. For the valleys we cal- 
culate the magnetic correlation for conductance, even in 
the crossover between ensembles of various symmetries. 

The QD attached to two leads, 1 and 2, is described 
by the Hamiltonian H — Hl + He, + iJLD, where 
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a=l,2 



corresponds to the leads, with vp being the Fermi veloc- 
ity, and cr =1, J. labels the spin. We consider only single 
channel contacts. The dot Hamiltonian is 
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where the first term describes the non-interacting dy- 
namics of the closed dot, and the second term, with 



corresponds to charging energy. Di- 
mensionless parameter A/" is a linear function of the gate 
voltage. The matrix H belongs to an RMT ensemble. 
In the case of quantum dots based on a 2D electron gas, 
the spin-orbit interaction is atypical and is characterized 
by two parameters and ey° . In the presence of a 
magnetic field with both perpendicular and parallel com- 
ponents to the gas plane (let the third axis be normal to 
the plane), Ti, can be written as 

^ \ Hoi + iX{xl + aj_a^) + ia{Aia^ + A2cr^) 



b^BhCT^] - ezl- a/2. 
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where Ho and Bh are real symmetric M x M matrices, 
with (TrHg) = M^ {TtBD = M2, and X and A are 
real antisymmetric matrices with {TiXX"^) = and 



(the limit M 

,2 



oo is taken event u- 



(TrA-4j) 

ally). Here = tte^/A, = 7re||°/A, is the Zee- 
man splitting energy due to the parallel magnetic field in 
the direction given by the unit vector I = {li,l2,0), and 
b\ — ne^/A where describes the combined effect of 
Zeeman splitting and spin-orbit scattering. The orbital 
effect of the magnetic field is characterized by the energy 
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e_B = x^A/tt = KETh{^/^o)'^, where k is a coefficient de- 
pendent on the shape of QD, $ is the flux of the magnetic 
field through the dot, and $o is the flux quantum. 

The tunnehng Hamihonian Hj^d couples the states of 
the leads a — 1, 2 to orbital states in the QD: 
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Santnci{k)c„ + h.C, 
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Here v — 1/{2itvf) is the leads density of states per spin. 

First we consider the statistics of the peaks. For the 
temperature range T <C A -C i^c, only the last occu- 
pied, possibly degenerate level participates in the elec- 
tron transport close to the peaks. Using the Golden 
Rule, the escape rates of the level a into the first and 
second leads are = -^91,2, where gi = A\ti\^yi is the 
dimensionless conductance of the z-th contact. Here 
9t 
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is a fluctuating quantity, and 4'a{i) are spinors with com- 
ponents (j,tJ. |Q!)- For the peaks we will confine our- 
selves to strong spin-orbit scattering ^ A. In 
this case, and for zero magnetic field, the random matrix 
Ti belongs to a Gaussian symplectic ensemble (GSE): it 
is a Hermitian matrix with quaternionic entries. In the 
notation of 11] this corresponds to the symmetry class 
(/3 = 4, E = l,s — 2). The number s — 2 signifies 
Kramer's degeneracy of energy levels. In the presence of 
a strong enough magnetic field such that either es 3> A 
or ez ^ A, the Hamiltonian Ti belongs to a GUE cor- 
responding to (/3 = 2,1] = 2,s — 1), and the symmetry 
group of the RMT ensemble is U(2Af). In this case, the 
levels are not Kramer's degenerate. In both cases, the 
random variables yi in Eq. (jSJ are statistically indepen- 
dent, and their distribution is given by 
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Assuming A 3> T ^ F" = {gi+g2)A/2Trh, one can use 
the rate equations In the GSE case, we will label 

the Kramer's degenerate levels by 4'a.s, with the s — ± 
states related by the operation of time reversal: V-'q- — 
Due to this relation, (V'l^+^/'a,+) = {ipi_tpa.-), 
and the escape rates are therefore equal for the s — ± 
states (see Eq. (O). In the absence of interactions, the 
level a can be in four states; empty, doubly occupied, or 
singly occupied with s — i. The picture changes when 
we take the interaction into account: due to the large 
charging energy, only three states with the number of 
electrons on the dot differing by 1 can be in resonance. 
For the peaks corresponding to M* = 2j + 1/2, (j G Z), 
the doubly occupied state has an extra energy Ec, and 
does not participate in transport (the case JV* = 2j — 
1/2 gives the same result for the peak height). In the 
stationary state the rate equations yield 
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A gig2 dfp/dx 
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where fpix) = 1/(1 + expx), and x = 2Ec{M - M*)/T. 
At the maximum which slightly deviates from x — Q 
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The only difference in the {f3 — 2,Y, — 2, s — 1) case is 
that now there is no Kramers' degeneracy and the reso- 
nant level can only be in two states: empty or occupied. 
Similarly to Eq. |(7J) we find 
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2 T gi+ g2 dx 
The maximum occurs at x = 0, and we find 
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Since the yis have the same distribution in both ensem- 
bles we see that Gpeak has the same distribution in both 
cases except for a scaling. Using the probability distri- 
bution Eq. ©, we obtain Gpeak in terms of the average 
contudctances ((71,2} and a single random variable a, 
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where Xi3=4,s=2 = 3 - 2^/^, and Xi3=2,s=i = 1/8. The 
probability distribution for a is given by (see Fig. ^ 
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Here Kq{x) and Ki{x) are MacDonald functions and 
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For the average we obtain 
2(1 + 62) 



1 - a' 



{a) = 8(1 - a)-! 
1 



105 



21 



-i^(6,3/2,ll/2;a)+ ^ 



F(5,l/2,ll/2;a) 
2 1/2 



(13) 



(14) 
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with F{. . . ) being the hypergeometric functions. Equa- 
tion ffT^ is well approximated by (a) ~ (8 — 3a)/10. 

Equation (|ll|l is our main result for the statisitical dis- 
tribution of the peak heights. We see that the application 
of the magnetic field causes the average conductance to 
drop by a factor 8(3 — 2'^/2) ~ 1.37, similar to antilocal- 
ization for bulk systems. Surprisingly, the shape of the 
distribution remains the same. This is because at strong 
magnetic field the RMT ensemble crosses over to a uni- 
tary ensemble with two channel leads corresponding to 
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FIG. 1: The probability distribution W{a) in different en- 
sembles, for symmetric contacts (a = 0). The f3 — 1,2, were 
studied in Refs. 



spin projections. The statistics of the eigenstates in the 
latter are the same as that of the symplectic ensemble, 
and hence the same distribution of peak heights. Though 
the drop in the average peak heights is a manifestation of 
the lifting of the Kramer's degenaracy, the numerical fac- 
tor for this drop is non-trivial. Had we ignored the charg- 
ing energy, we would instead obtain a reduction by a fac- 
tor of 2. The different numerical factor originates in the 
exclusion of the doubly occupied degenerate level from 
transport by the electron-electron interaction. We em- 
phasize, however, that the distribution function changes 
for intermediate values of the perpendicular field. 

Next we turn to the mesoscopic fluctuations of the val- 
leys, \Af -J\f*\ > T/Ea. Then, Eqs. 0-® suggest an 
exponentially small conductance. However, this is incor- 
rect, as the rate equations are based on processes of the 
first order in Hld. For such processes a charge is trans- 
fered between a lead and QD, and the conductance in the 
valleys is small because of the charging gap in the final 
state. This is not the case for higher order processes that 
allow for both the initial and the final excitations to be 
in the leads, so there is no gap for the final state. Such 
processes are referred to as co-tunneling 0, , and are 
suppressed only algebraically, G{N) oc l/{2Ec\N -N*\). 
We consider elastic co-tunneling which is the dominant 
process in the temperature range T <C (AE'c)^^^- In 
these processes a charge is transferred from lead to lead 
via a virtual transition to an excited state in the dot. The 
conductance calculated using the Golden rule is 0, Q 
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where the amplitudes A'^^"''^ i^V''^^) correspond to pro- 
cesses in which an electron(hole) in spin state cti in the 
first lead is transferred to the second lead in spin state 
a2- The second order perturbation theory in Hljj gives 
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where a is summed over energy levels of the QD. Here 

Ee = 2Ec{N* - N), Eh = 2Ec{Af - Af* + 1), (17) 

(TV* - 1 < JV < A/"*), are the electrostatic part of the 
energy of the virtual state for the electron-like and hole- 
like processes respectively. The eigenenergies are mea- 
sured from the last occupied level so that the step func- 
tions 6{±ea) select empty(filled) states for electron(hole) 
Hke processes. Equation l|15|l can be expressed in terms 
of the Green functions (GF) for the closed dot: 
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de gi^2,a-itT2(g) ~ ^i2,<Tia2(e) 



Using this relation we can express all the moments of 
the conductance in terms of the average of a product of 
GFs. Furthermore the condition A <^ Ee,Eh, allows us 
to use the diagrammatic technique in terms of diffusons 
and cooperons to calculate the latter. In this approxi- 
mation the GFs become Gaussain variables with average 
{Qnrif^) <x. 5nm, and variances givcn in terms of the Dif- 
fuson and Cooperon matrices 0, Q 
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Tr (7^gf2(ei;Si)a''gf2(e2;i?2)) = 
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where lo — ei 



£2: 



cr''"*' = i, and cr', i = 1,2,3 are 



Pauli matrices. Here we defined ^/l = (^^1^2) o-^- For 
the Hamiltonian Eq. JSJ, the inverse of the diffuson and 
cooperon matrices are given by [TH Il2l | 
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C-\uj; Bi,B2) = -iuji + it^f] + 
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(21) 



where ff^ — —ie^^^, and 77^^ — lf_iSi^,o + ^iy<5^,o, with / the 
unit vector in the direction of parallel field, and 
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Bj + B^ 
2$oM 



2$oM 
where A is the dot's area. 

According to Eq. ifT^ll the amplitudes Fe,Fh are lin- 
ear in and and are also Gaussian variables with 
zero average, so that we can calculate the moments of 
the conductance, Eq. (fl^ . using Wick's theorem. From 
Eqs. (fH^ and (pn|l for ^ = = 0, we obtain the result 
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FIG. 2: The correlation function C12 for By = 0, = 0.1, 

and ax = ^J^fjE = 0, 1, 2. The function is non-analytic at 
AB/Bc = 0, ±2a, where Bc^l = ^oV-B/K^Th. 



independent of all the crossover parameters. Using 
Eqs. (1201)- (ED together with ((TH|l . we calculate the corre- 
lation of the conductances at different values of the per- 
pendicular magnetic field. In the vicinity of the peaks 
corresponding to the condition Eh -C or i?e <C Eh, 
for C12 = {5G(Bi)5G[B2)) / {G)^ we obtain 



Furthermore for strong spin-orbit scattering, e^, ej*" ^ 
E, we can calculate all the moments of the conductance 
G, and obtain the distribution functions for ^ — Gj (G): 
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, (26) 
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A^), <f> is the magnetic 



where A = A (<i>2/$^), 7 = 7/(1 
fiux, and <i>c = <i>o ^/Ejn^Th- 

In conclusion, we studied the statistics of the conduc- 
tance for CB peaks and valleys in the presence of spin- 
orbit scattering. We calculated the distribution function 
of the peak heights for strong spin-orbit scattering in 
the presence and absence of time-reversal symmetry. We 
found that the average peak height is reduced in the latter 
case. For the valleys we calculated the average conduc- 
tance and the correlation function of conductance for a 
2DEG QD, as a function of perpendicular magnetic field. 

We are grateful to C. M. Marcus for suggesting this 
problem. 
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where A(z) = -^{z) — A(— z) and 
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for complex z, with Li2(x) being the dilogarithm func- 
tion. The asymptotic behavior is A(z) = 1 + zlnz/ir, for 
\z\ <C 1, and A(z) — {Trz)~^\n^z, for Rez 3> 1. Here 
E ~ min{Ee, Efi), and A„ and A„ are the four eigen- 
values of 2?~^(0; Bi, B2) and C~^(0; B2) matrices re- 
spectively. The 4x4 matrices Aofi are given by 



Ad = U^^LdUd, Ac = Uc^LcUc, 



(25) 



where Ud and Uc are the matrices whose columns are the 
eigenvectors of 25^^(0; Bi, ^2) and C~^(0; Bi, ^2) respec- 
tively, L^^j^ = 6t"'-2S^''S^'' andL^" = S>"' -26°t'6°'' . For 
zero parallel field and small perpendicular fields Zeeman 
splitting can be ignored and the eigenvalues of inverse 
cooperon and diffuson correspond to the usual singlet 
and triplet states, and we have (see Fig. 
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